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ABSTRACT 

A simplified one-dimensional model of reactive flow is 
presented which captures features of aeropropulsion sys- 
tems, including acoustically driven combustion instabil- 
ities. Although the resulting partial differential equa- 
tions are one-dimensional, they qualitatively describe 
observed phenomena, including resonant frequencies and 
the admission of both steady and unsteady behavior. 
A number of simulations are shown which exhibit both 
steady and unsteady behavior, including flame migration 
and thermo-acoustic instabilities. Finally, we present ex- 
amples of unsteady flow resulting from fuel modulation. 
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INTRODUCTION 

The current drive toward lean operating, low-emission, 
low-loss combustors in next generation aeropropulsion 
systems makes them susceptible to various instabilities 
and makes the understanding of those instabilities crit- 
ical. Combustion instabilities, specifically acoustically 
driven instabilities, often arise from the interaction be- 
tween the energy released by the chemical reactions and 
the underlying fluid flow. This coupling can dramatically 
alter the dynamical behavior of the combustion process, 
and result in unstable behavior. Typical combustion in- 
stabilities in propulsion systems can often be character- 
ized by oscillatory pressure fluctuations coupled with an 
unsteady heat release. This instability can lead to flame 
blowout, reduced performance, and decreased engine life. 


1 



Although the magnitude of the pressure oscillations 
can be substantial, the total energy associated with the 
unsteady behavior is typically only a small fraction of 
the total energy generated by the combustion process. 
As a result, the pressure fluctuations can often be iden- 
tified with the acoustic modes of the enclosed volume. 
The acoustic behavior of the propulsion system induces 
an oscillatory release of heat, which in turn drives the 
acoustic modes, or vice versa. 

The fundamental equations of reactive gas dynam- 
ics can be described as a set of reaction-diffusion equa- 
tions governing the energy release associated with the 
combustion and transport of the chemical species, cou- 
pled to the Navier-Stokes equations for the underlying 
fluid flow, subject to an equation of state. Because of the 
complexity of both the fluid and transport equations, as 
well as the coupling mechanisms, the ability to model 
and accurately predict the dynamical behavior and per- 
formance of propulsion systems is often limited. As a 
result, the simulation and characterization of unstable 
combustion behavior is critical to the development of 
future propulsion systems and is the subject of intense 
interest [1]. 

A common approach to modeling thermo- acoustic 
combustion instabilities is to simplify the governing 
equations to those of one-dimensional flow [1, 2]. The 
use of a constitutive model with a sensitive time-lag for- 
mulation is then employed to couple the unsteady heat 
release with the underlying fluid flow [3] As a result, 
the time-lag, or n-r formulation, eliminates the trans- 
port equations which govern the chemical kinetics, and 
reduces the complexity of the resulting model. Although 
this reduction presents a simplification of the model, it 
has met with great success in characterizing these insta- 
bilities. Despite this success, we study these instabilities 
with the inclusion of simplified chemical kinetics. Our 
goal is to obtain fundamental insight into the behavior 
of these systems, in particular the coupling between the 
unsteady heat release and the acoustic response. 

This work presents a numerical investigation of 
both the convection and transport equations governing 
one-dimensional reactive flow. These equations repre- 
sent a simplified model for the flow in aeropropulsion 
combustors, and hence cannot capture such phenom- 
ena as flow oscillations driven by shear layer instabil- 
ities or the liquid-fuel vaporization process. However, 
they nonetheless exhibit self-excited thermo-acoustic in- 
stabilities and can provide a tool for studying the quali- 
tative behavior of such systems, including the sensitivity 
of the system to parametric variations in the convection 
and transport parameters. In particular, we study the 
dynamics of this system with an Arrhenius model for the 
chemical kinetics and a description of the flame holder. 
The simulation will be described, including the governing 
equations, justification for the simplifications, and inte- 
gration techniques. The results will be presented and 
contrasted with respect to the development of instabili- 
ties. 


GOVERNING EQUATIONS 

We consider a premixed gas flowing in a finite length 
combustor with constant cross sectional area, which un- 
dergoes a single-step irreversible chemical reaction. In 
this formulation, we assume that the mixture contains 
a single deficient component, the amount of which con- 
trols the reaction. In addition, this reactive component 
is assumed to exist in an inert bath, in which the total 
concentration of the reactive component is small — as the 
combustion process occurs, the overall gas parameters of 
the mixture remain unchanged. 

In a single spatial dimension, the basic nondimen- 
sional equations of motion governing reactive gas dy- 
namics are: 


U t + F x = 5, 


(i) 


with: 
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E represents the total energy and is defined as: 


E - ^ +/i 2 ^ 7 - 


Time is n or di mensio n alized by the acoustic time scale, 
so that p = 1.00, and Z(.) = dz/d( ). These equations 
are completed by the ideal gas law, p = pT. 

We assume that there is no external heat transfer, 
although these equations will be augmented below to 
include the effect of a flame holder. These equations are 
written in conservation form, with p(x, t) and u(x, t) the 
total density and mass weighted velocity, respectively, of 
the gas. The temperature is described by T(x, t), p(x, t) 
is the press ire, and y(x, t) represents the mass fraction 
of the react mt in the combustible gas. 

The sc urce terms within 5 contain the thermody- 
namic effec s of the reactive flow, as well contributions 
from diffusional terms. Q represents the rate of heat 
release by "he chemical reaction, while m is the mass 
conversion "ate of the unreacted fuel. These equations 
have been used by a number of authors to study flame 
propagation through high activation energy asymptotics 
[4, 5, 6, 7, 8] and the behavior of premixed combustion 
in wave rotors [9]. 

The r ondimensionalization has been performed 
with respect to the unreacted gas and the dimensional 
combustor length /. At the inlet, x = 0, we impose 
t/( 0, t) = 1 The boundary conditions on the pressure 
and tempei ature are dependent on the geometry of the 
combustion chamber. In addition, physically we expect 
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y(x,t) -+ 0 as x L It should be noted that the diffu- 
sive time scale is not necessarily controlled by the fluid 
properties but may well be controlled by the flowfield 
(e.g. through turbulent diffusion). 

Heat-release 

We describe the instantaneous rate of heat release, Q , 
associated with the combustion process as a one-step, 
irreversible reaction, where Q is proportional to the mass 
conversion rate of the unreacted fuel. In addition, we 
assume that m is proportional to the mass of unreacted 
fuel, with a temperature dependent reaction rate: 

Q = 0m, 

= 0 ((Apy) • fc(T)) , 

where 0 is the nondimensional heat release parameter 
and m is the mass conversion rate of the unreacted fuel. 
p and y are the density of the gas and the mass fraction 
of the premixed fuel, respectively, so that the quantity 
py physically represents the mass per unit volume of fuel 
available to the combustion process. A is a parameter 
associated with the rate of mass conversion and k(T) is 
the temperature-dependent component of the reaction 
rate. 


With this choice for k(T), the ratio between the maxi- 
mum rate of reaction at T = T\ and the initial rate of 
reaction at T = Ti is approximately: 


fe(Ti) 

k(Ti) 



Tq + 0 


)}• 


Flame Holder 

In the absence of a flame holder, a stable flame requires 
a delicate balance between the upstream diffusion of the 
flame and the downstream convection of the fluid. As 
a result, for most ail gas parameters, this balance will 
not occur and the flame will either propagate upstream 
to the inlet (flashback), or propagate downstream and 
out of the combustor, a condition known as blow-out. 
The flame holder acts to limit the propagation of the 
flame upstream. We note that flow effects, such as a 
spatially varying flow velocity due to cross-sectional area 
variations, can also serve to anchor the flame. 

We assume that the effects of the flame holder are 
twofold: (i) the diffusion constant S becomes spatially 
dependent, resulting from increased turbulent diffusivity 
downstream of the flame holder, and (ii) the flame holder 
transfers heat to and from the fluid. We incorporate 
both effects in this model. 


Mixture thermodynamics. The thermodynamic proper- 
ties of the premixed gas are contained within the pa- 
rameters 0 and A. If we subsequently consider varying 
mixtures of oxidant and reactant, both parameters are 
limited above by some maximum values, say An&x and 
Anaax- 

The flame temperature resulting from the combus- 
tion process scales with 0, the nondimensional heat re- 
lease, so that the temperature at the exit is approxi- 
mately T\ = T(l, t) = To + 0, where T 0 = T(0,t). 
As such, for a given reactant, the maximum value of 
0 roughly corresponds to stoichiometric conditions, so 
as 0 approaches /3 raa x , the equivalence ratio of the mix- 
ture approaches unity. Also, A, the Damkohler number, 
scales the rate of the chemical reaction. These quantities 
parallel the equivalence ratio and the heat content of the 
premixed fuel. 


Chemical kinetics. In this simulation we model the 
chemical reaction by a single-step, first-order Arrhenius 
model, so that the temperature dependent reaction rate 

is: 


k(T) 


A u . _ e{i-i/T) 

Arrhenius — c i 


with 6 the activation energy [4]. In addition, we also 
incorporate an ignition temperature Ti, below which no 
combustion occurs. As a result, k(T) takes the form: 



0, 

e *U-l /T) 


T <Ti, 
T > Ti. 


Diffusions I effects. We assume that the flame holder 
acts to increase the diffusion due to induced turbulent 
mixing. This implies that the scale of the turbulent ed- 
dies is larger than the width of the laminar flame zone. 
As a result, we model the spatial distribution of <5(x) as: 

6{x) = Si ^ ^1 - (tanh(20(x - a)) + l) + ~ . 

In the limits x — ► ±oo the diffusional coefficient ap- 
proaches So and Si respectively. 


Thermodynamic effects. The heat transfer effects are as- 
sumed to be spatially dependent as well, of the form: 

/thermal (T(x, t), x) = 9{ x ) (T f ~ T(x ,t)^ , 

= / o exp[(l-(50(x-a)) 3 4 )] X 

( T/-T(*,0), 

where g(x) is the local convection transfer coefficient, 
and Tj is the nominal temperature for the flame holder 
[10]. Thus g{x) describes the spatial distribution of the 
convective effects of the flame holder. We augment the 
energy equation, by /(T,x). In addition to the increased 
diffusional effects discussed above, the flame holder acts 
to limit the propagation of thermal energy upstream of 
the flame. This form for the heat transfer coefficient 
provides an analytic approximation which is localized in 
the neighborhood of x = a. 
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NUMERICAL SIMULATIONS 
With prescribed initial conditions, the system governing 
reactive gas dynamics in the form of Eq. (1) is simu- 
lated using MacCormack’s method [11], which is second- 
order accurate in both space and time. The flow is 
solved in the interval 0 < x < 1 for t > 0, subject to 
characteristic-based boundary conditions [12], and the 
diffusional terms are incorporated as a source using a 
second order spatial differencing scheme. To verify the 
qualitative accuracy of our integration scheme we have 
varied both the number of grid points, as well as the size 
of the time step. 

Unless noted, we assume 7 = 1.40 and choose the 
flow parameters: 

Sc = 1.00, Le = 1.00 , /a = 1.00, 

0 = 1.50, 6 = 2.00, Ti = 1.00. 

The flame holder located to be at 25% of the combustor 
length (o = 0.25), with: 

fo = 0.1, Tf = 1.00, Si = 0.01, ^ = 0.025. 

Ol 

In the following, we present several simulations to 
illustrate the dynamical characteristics of the model, in- 
cluding both steady and unsteady behavior. In addition, 
we consider both choked and open inlet conditions, as 
well as a number of different parametric conditions. In 
the figures, darker shades indicate lower amplitudes. 

Choked-open Boundary Conditions 

The inlet conditions are choked and isothermal, so that: 

p( 0, t)u( 0, t) = 0.20, T(0, t) = 0.75, 

while the exit is open: 

p(l,t) = 1.00. 

The remaining boundary conditions are determined from 
the interior of the domain of integration x € [0, 1]. In 
Figure 1, we show the steady response for A — 2.50. We 
note that this response is stable. 

Flame Migration 

By starting the flame away from its equilibrium posi- 
tion, we can also observe the transient behavior of the 
flame. As seen in Figure 2, combustion initially occurs 
at x = 0.40, well away from the flame holder located 
at x = 0.20. The simulation parameters are chosen as 
above, with A = 1.50. As time evolves, the flame moves 
upstream, eventually stabilizing on the flame holder. 

Open-open Boundary Conditions 

We now apply open boundary conditions and assume the 
flow is isothermal at the inlet, so that at x = 0: 

p( 0, t) = po = 1.05, T(0, t) = 0.90, 

while the exit is open: 

P(M) = 1-00, 


which implies a 5.0% pressure drop across the combus- 
tor length. The remaining boundary conditions are de- 
termined from the interior of the domain of integration 
*€[ 0 , 1 ]. 

Unsteady Response 

For A = 1.50, the asymptotic solution is unsteady, and 
the fundamental mode is excited (see Figure 3). In Fig- 
ure 3a we show the time trace of the pressure at x = 0.25. 
As time increases the amplitude of the pressure oscilla- 
tions initially grow before limiting around t = 35. For 
comparison with Figure 5, in Figures 3c, d we show the 
steady solution, which for this value of A is unstable. 

In Figure 4 we show the time-varying components 
of the flow, denoted by (•)', including the unsteady rate 
of heat release <j'(x,t). We define the local Rayleigh 
index 7£(x,t) = Q'(x,t) -p'(x,t), and further define the 
Rayleigh index as the time average of this local index 
over one oscillation: 

G(x)=^ 

where T is the period of oscillation. 7 Z(x,t) represents 
the work done to the fluid by the chemical heat release at 
time t and spatial location x. As time evolves the local 
Rayleigh index, and hence the work done, takes on both 
positive and negative values. However, G(x) represents 
the average work done over one cycle. The well-known 
Rayleigh’s criterion states that if G(x) > 0 energy is 
added to the system and the amplitude of the pressure 
oscillation is expected to increase, while for G(x) < 0 
the flame-acoustic interaction removes energy from the 
system. As expected, examination of Figure 4 indicates 
that over most of the flow field the unsteady heat release 
is in-phase with the oscillating pressure component — the 
Rayleigh index is positive and we indeed find an un- 
steady response. 

Parametric Variations 

Darling, et al. experimentally studied the acoustic re- 
sponse in several combustor configurations [13]. In then- 
work, it wa. noted that "... the resonances were very 
sensitive to changes in the mean flow conditions, and the 
effects of these changes were not always consistent.” For 
example, at low inlet temperatures the amplitude of a 
resonant acoustic mode was increased as the equivalence 
ratio was increased. However, for larger inlet tempera- 
tures the resonant peak was decreased as the equivalence 
ratio was in :reased. 

Although, as mentioned above, we do not attempt 
to make diiect comparisons between our fuel parame- 
ters (0, A) ;md the equivalence ratio and heat capacity, 
we do note that as we change 0 and A this model pos- 
sesses similar qualitative behavior. For both A = 0.75 
and A = 4 00 the asymptotic response, shown in Fig- 
ure 5, is steady, while, as seen in Figures 3 and 4, an 
unsteady response exists for A = 1.50. This type of 
behavior is common — there is an intermediate range of 
A for which the steady-state solution is unstable. For 
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a) Utq(x)) vs. x b) (T c ,(x),y eg (x)) vs. x 

FIGURE 1: Choked-open boundary conditions — steady response (A = 2.50). The flow velocity and temperature are 
shown dashed in a) and b) respectively, while the pressure and species fraction are solid curves. 



a) 


b) 


FIGURE 2: Flame migration and stabilization with choked-open boundary conditions. Figure b) is an expanded view 
for 0 < t < 2 which emphasizes the migration of the flame, initially located at x = 0.4, upstream to the flame holder 

at x = 0.2. 
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a) p(0.25, t) vs. t 


b) (p(0.25, t), X(0.25, t)) vs. t 



C) (p e<? (x),Ue«(x)) VS. X 


d) (^(^),y«g(x)) VS. X 


FIGURE 3: Open-open boundary conditions (<5i = 0.02). For A = 1.25 tae response is unsteady. Figure a) presents 
the growth of the pressure response as time increases, eventually approaching a periodic oscillation. In b) the steady- 
state oscillation in pressure and temperature at x = 0.25 are shown as the solid and dotted line respectively, while 
in c) and d) we show the time-independent profile which is unstable. The flow velocity and temperature are shown 
dashed in c) and d) respectively, while the pressure and species fraction ire solid curves. 
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c) u'(x,t) 


d) 


FIGURE 4: Open-open boundary conditions (<5i = 0.02). Unsteady components of the response shown in Figure 3. 
The time interval shown is while the amplitudes of the oscillations are small. As seen previously in Figure 3, large 
amplitude oscillations eventually develop. 
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A sufficiently small or large, the steady-state is stable. 
Reflecting back on our discussion between the relation- 
ship between 0 and A, and the equivalence ratio, this 
implies that for near-stoichiometric mixtures, our model 
leads to an unsteady asymptotic response. As A -4 0 
we numerically encounter flame blowout, while large A 
admits the possibility of flashback. 

Finally, we note that for 0 = 2.0 the ratio between 
the final and initial burning rate is approximately 3.2. 
If we decrease the activation energy to 0 = 1.0, so that 
ratio between the final and initial burning rate is ap- 
proximately 1.8 the response is stable for all values of 
A. 

Fuel Flow Modulation 

Looking forward to possible attempts to implement con- 
trol strategies based on fuel modulation, we study the 
effects of small amplitude oscillations in the fuel species 
parameter y at the inlet: 

y(0, t) = 1.0 + 0.01 sin(27ru>t). 

In Figure 6a, we show the amplitude of the pressure 
oscillation at x = 0.25 as the modulation frequency u 
varies from u; = 0.25 to u = 1.00. When the fuel species 
is modulated near the fundamental acoustic frequency 
at w ~ 0.6, modulation amplitudes of 1% lead to large 
amplitude acoustic oscillations. Figure 6b presents the 
growth of the pressure oscillation for u = 0.6. We note 
that in these simulations the flame zone is located at ap- 
proximately 20% of the combustor length. As a result, 
diffusional effects reduce the amplitude of the fuel mod- 
ulation at the flame zone when compared to the inlet 
fluctuations. Thus we expect that if the fuel was modu- 
lated slightly upstream of the flame zone, rather than at 
the inlet, the resulting acoustic response would increase 
further. 

SUMMARY 

This work focuses on numerical simulations of reactive 
flow which, while retaining the simplification of one- 
dimensional flow, contains equations for species trans- 
port as well as fluid convection and diffusion. As evi- 
denced by the representative figures shown, it is possible 
to study thermo-acoustic instabilities, flame evolution, 
and the behavior resulting from parametric variations. 
In addition, the model for the chemical kinetics contains 
two effects — Arrhenius temperature dependence and an 
extinction temperature. We also investigate the dynam- 
ical behavior with the inclusion of both thermodynamic 
and diffusional effects from the flame holder. 

Examination of the unsteady temperature distribu- 
tion indicates two distinct regions along the width of the 
combustor — i) the flame zone, in which large tempera- 
ture variations occur due to the chemical reactions, and 
ii) remainder of the combustor, in which the temperature 
variations are, in fact, almost isentropic. 

Although the model is too simplified to obtain 
quantitative comparisons with the phenomena observed 
in experimental tests, it provides an ideal framework to 


study the qualitative behavior of thermo- acoustic cou- 
pling and can be used as a basis for the investigation of 
the control of these oscillations and the implementation 
of various control strategies, including modulation of the 
fuel flow. 
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a) A = 0.75 


b) A = 4.00 


FIGURE 5: Open-open boundary conditions — for A = 0.75 and A = 4.00 the time-independent response is stable, 
shown in a) and b) respectively. As A increases the length of the flame zone decreases. In each figure the states axe 
denoted by: pressure ; temperature ; velocity ; species 




a) i/(0, t) = 1.0 + 0.01 sin(27ru;t) 


b) p(0.25, i); u) = 0.60 


FIGURE 6: Acoustic oscillations resulting from fuel flow modulation (0 = 1.50, x ~ 0.25). In a) the amplitude of 
the response is plotted versus w the forcing frequency, while in b) the pressure is shown for = 0.6 — the initial 
conditions in b) coincide with the steady response for no fuel modulation. 
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